Scaling behavior in Spiral Defect Chaos. 
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We find the evolution toward power-law scaling in the distribution of roll lengths and nearest- 
neighbor distributions in a weakly turbulent regime of Rayleigh-Benard convection, known as spiral 
defect chaos. The state has a bounded domain of wave vectors in Fourier space attributed to the flow 
being highly confined to two dimensions. Our results indicates the existence of power-law scaling in 
the unconstrained horizontal dimension. The techniques described are broadly applicable to other 
pattern forming systems as well. 
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The patterns exhibited by spatially extended systems 
reflect the underlying dynamics. In non-equilibrium sys- 
tems, complex patterns emerge at a macroscopic scale 
as a result of local nonlinearities Relating the dy- 
namics from the microscopic scale of nonlinearities to the 
macroscopic patterns observed remains challenging. It is 
therefore important to develop generally applicable char- 
acterizations of a system that can capture the dynamics 
across a multiple scales in a robust and consistent way. 

The emergence of patterns represents the delineation 
of the system into regions with different physical charac- 
teristics. For instance, the domain structure of a mate- 
rial distinguishes between regions of differing orientation. 
The different regions may be considered as the different 
components comprising the pattern. In fluid systems, 
these structures are often in the form of variations in 
local velocities during the flow. For example, a charac- 
teristic of turbulent flows is the formation of eddies at a 
multitude of scales, in contrast, laminar flows have a well 
defined length scale associated with them 0, 0] • 

The flow of a fluid in highly constrained geometries is 
dominated by the confining boundaries. The reduction 
in size of a system often results in the suppression of tur- 
bulence through an effective reduction of the Reynolds 
number associated with the flow [3(. While this reduc- 
tion in size can be imposed along a single spatial dimen- 
sion without affecting the other spatial dimensions, the 
influence on the resulting flow is not easily extracted. 

In this letter, we extract power-law scalings associ- 
ated with patterns representing fluid flow in Spiral De- 
fect Chaos, a weakly turbulent state of Rayleigh-Benard 
convection 0, [B|, 0, 0- In this flow, the vertical di- 
mension of the flow is highly constrained in comparison 
to the horizontal dimension. The convective instability 
is coincident with the formation of roll like structures 
with a wavelength on the order of the depth of the con- 
vective cell [H, The system has received much at- 
tention as it shows unexpected chaotic behavior that is 
bistable with at values of system parameters correspond 
ing to the flow being stable and stationary j^, flQ , 




FIG. 1: A typical image of the convective state exhibiting 
spiral defect chaos(at e = 2.63) as visualized using shadowg- 
raphy. The convective fluid is bounded horizontally in a circu- 
lar domain with radius 3.8cm and vertically extends to about 
0.69mm 



tration and mean flows that may be linked with defect 



nucleation [H, E, Q [H, E 



A number of characterizations of the state have been 
made emphasising local instabilities, wave-vector frus- 



We use two characterizations of the state that reveal 
length scales that are associated with the structure of 
the flow in spiral defect chaos. The important differ- 
ence with earlier characterizations is that these measures 
reveal underlying power-law scalings in real space not 
identifiable with the distributions of the Fourier power 
spectrum. In the first characterization, we determine the 
scaling of the area of distinct convection roll structures. 
In the second characterization, we reduce the patterns 
exhibited by the system to a graph. This parameterizes 
the system through number of nearest-neighbors of each 
distinct roll. 

Our experiments use compressed carbon di-oxide as 
the convective fluid. The gas is held in a convection cell 
bounded below by a 1cm thick gold plated aluminium 
mirror and above by a 2cm thick transparent sapphire 
crystal. The lateral boundaries comprise of stacks of fil- 
ter paper 0.069cm thick that also act as a spacer between 



the upper and lower boundaries. The filter paper has a 
3.8cm diameter circular hole in the center within which 
the convection is visualized. The aluminium mirror has a 
resistive thin film heater attached to it on the bottom and 
rests on high precision screws that may be turned to align 
the mirror. The mirror is levelled with respect to the sap- 
phire window using laser interferometry. Chilled water at 
a controlled temperature is circulated over the top of the 
sapphire window to maintain its temperature. In addi- 
tion, thermistors embedded in the side of the aluminium 
mirror as well as close to the sapphire plate are used to 
monitor the temperature gradient across the convective 
layer. The setup in contained within an aluminium canis- 
ter pressurized to 30atm with CO2 gas using an external 
pressurized cylinder. Computer control of the thin film 
heater is used to set the temperature difference, T, be- 
tween the top and bottom boundary of the convective 
cell to within lOmK of the desired value. Our experi- 
mental setup is similar to that used earlier by Rogers et. 
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At a critical temperature difference across the fluid, 
T = T c , the onset of convection occurs resulting in the 
formation of roll like structures. These rolls represent the 
transport of fluid vertically across the cell when thermal 
expansion induced buoyancy in the fluid overcomes dis- 
sipation due to viscosity and thermal conduction [3], Q . 
The visualization of the convective rolls is done using 
shadowgraphy 1|J 2^] ■ This technique utilizes the varia- 
tion in the refractive index between hot and cold regions 
of the gas to focus light differently. Light from a bright 
halogen lamp is passed through a pinhole and collimated 
using a concave lens. After passing through the convec- 
tive layer, the spatial variation in light intensity is cap- 
tured using a CCD array. The image is indicative of the 
changes in the refractive index of the light as it passes 
through the fluid. In our setup, bright and dark regions 
represent hotter upflows and cooler downflows. 

We set our experiment at different values of e(= (T — 
T c )/T c ), and hold the temperature for over one hour to 
let transients dynamics die out. A sequence of 8000 im- 
ages acquired at 11Hz are used in our data analysis. The 
total time the dynamics of the system is captured is many 
times the vertical thermal diffusion time (about 2 seconds 
in our system) . The response of system falls qualitatively 
into three regimes [U, Ejj: low, < e < 0.5, where the 
convective rolls are stationary states of the system, in- 
termediate values of 0.5 < e < 1.5 result in the onset of 
time dependence and weak turbulence, and at 1.5 < e, 
the system shows prominent time-dependent behavior. 

A background image of the fluid flow prior to the on- 
set of convection is subtracted from every image to re- 
move stationary optical inhomogeneities that may be as- 
sociated with irregularities in the optical path(such as 
scratches in the mirror etc.). These difference images 
are Fourier filtered to eliminate the influence of high fre- 
quency noise associated with the CCD array as well as 




FIG. 2: In figure (a), a set of convection rolls representing 
the dark areas of figure [l] are indicated by dark regions of a 
binary image. As an example of a distinct component of this 
flow, one of the rolls that is singly connected is shaded in gray. 
In figure (b), the complementary flow is shown, representing 
the bright areas in figure [T] In addition, the gray component 
is added for visual effect, while also indicating its nearest 
neighbors in color. 



low frequency variations associated with nonuniform il- 
lumination. The intensity profile of the resultant image 
is representative of the structure of the convective rolls. 
The image is thresholded at the median value of intensity 
to yield a binary image. This choice of a threshold segre- 
gates the convective region into equal areas representing 
hot upflows and cold downflows to be consistent with the 
conservation of mass. 

The segregation of the flow into two distinct compo- 
nents as depicted by the binary image (figure [2]) provides 
the basis for our characterizations. As illustrated in fig- 
ure each distinct bright (dark) roll is completely sur- 
rounded by dark(bright) rolls, or the boundary of the 
convective cell. In the first characterization, we simply 
look at the area of distinct bright or dark rolls as a frac- 
tion of the total convective area. A histogram of the 
number of distinct rolls at a given area fraction divided 
by the total number of distinct rolls vs their relative area 
at a given value of e is plotted in figure [3] Since the con- 
vective rolls have a well defined average width, the area of 
the rolls is well approximated by a constant scalar multi- 
ple of the length of the roll. The convective rolls may be 
looked upon as primarily one dimensional curves within 
the circular convective region. 

When e > 1.5, a power-law scaling of the roll areas to 
number of rolls emerges. However, prior to the occur- 
rence of this scaling, we note two distinct scalings occur 
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FIG. 3: The area of convective rolls scale as a power-law with 
the number of rolls. The slope in the plot (a) for both upflows 
and downflows is close to —1.8. For comparison, in plot(b), 
the curves are split into regions of slopes approximately —0.35 
and —3.6 at e — 0.87. Upward pointing triangles refer to 
bright regions, and downward pointing triangles to dark re- 
gions of the binary images captured. 



based on the area fraction of the rolls as seen in figure [3] 
The kink in these plots occur at an area corresponding 
to that of a roll with length equal to the diameter of the 
convective region. There is a sharp drop in the number 
of rolls that have a length greater than the diameter of 
the convection cell. This is because the rolls are without 
much curvature at lower values of e. At higher values of e, 
the roll curvature increases, with an associated increase 
in the number of rolls with a larger area. The distinction 
between the distributions for upflows and downflows sug- 
gests that the flow is non-Boussinesq in nature fH I22I [23j| . 
The asymmetry between upflows and downflows increases 
with an increase in e. We find the system to be composed 
of a higher number of dark rather than bright rolls. 

A further topological characterization may be made 
by considering the binary image of the flow as a bipar- 
tite graph. In order to do this, we consider the number 
of distinct bright(dark) rolls each distinct dark(bright) 
roll shares an interface with. Figure EJd illustrates this 
measurement by indicating the nearest neighbors of a 
single convective roll in the flow(indicated in gray). In 
our binary image, bright (dark) regions are surrounded 
by dark(bright) regions with the only exceptions being at 
the boundaries. We count of the number of nearest neigh- 
bors every dark(bright) region has that is bright (dark). A 
set of numbers representing the number of nearest neigh- 
bors for each distinct dark/bright region is thus gener- 
ated at a constant value of e. We normalize this set so 
their sum is unity to compute the probability of hav- 
ing a given number of nearest neighbors. This is ploted 
against the number of nearest neighbors for a convec- 



tive roll as shown in figure [4j The probability of having 
a given number of neighbors scales as a power-law with 
the number of nearest neighbors. We see no significant 
difference between upflows and downflows in the scaling 
of the number of nearest-neighbors. This measurement 
is therefore insensitive to the breaking of the Boussinesq 
symmetry in contrast to the results presented in figure [3] 
Figure 2] shows the scalings for bright regions of figure [2ji 
and the results for the dark regions do not show signifi- 
cant differences. 
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FIG. 4: The scaling of the number of nearest neighbors with 
increasing temperature gradient across the system. At high 
values of e, the system evolves to show a power-law scaling. 
The slope at the highest value of e shown above is —2.9 



The probability of having only one neighbor is finite for 
rolls that are at the boundary, or are completely encircled 
by a different roll. We note the general trend in figure [4] 
for the number of rolls with two nearest-neighbors to de- 
crease with increasing e. In the simple geometries of par- 
allel rolls and concentric circles, the number of nearest- 
neighbors is two. In single spiral formations, the number 
of nearest neighbors would be one, there being only a sin- 
gle bright roll and a single dark roll. When more complex 
patterns emerge, they can often be broken down into spa- 
tially localized regions with simple geometries [2j|. The 
present characterization develops a general technique to 
make global characterizations of the system based on the 
nearest-neighbor distributions within the flow. 

The experimental characterization we have used relies 
on a small image series at each value of e. It has been 
shown that the state of spiral defect chaos exhibits tran- 
sitions between straight-roll states and curved roll states 
that are time dependent. This was done by observing 
the evolution of the system over a period of a month, 
keeping e fixed Uj . The importance of our study lies in 
extracting global measures that reveal an underlying scal- 
ing behavior for the various patterns exhibited. A more 
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accurate parameterization of this scaling with e would 
involve averaging over very large experimental data sets. 
Our current data suffices to demonstrate the existence of 
power-law behavior. 

While the lengths of rolls show a sharp kink that is 
associated with the system size, there is no feature cor- 
responding to this in the nearest-neighbor distributions. 
The scalings of the area of the system are metric de- 
pendent, and correspondingly show an associated length 
scale of the system size. The nearest-neighbor distribu- 
tions however are metric independent and a purely topo- 
logical measure. The two measures are related as they 
determine the packing of upflows and downflows in the 
system. A larger packing fraction of the rolls translates 
to an increase in the number of rolls. The increase in the 
number of distinct rolls with an increase in e has been 



reported earlier [22j, |23[ . 

The techniques presented allow for connecting dynam- 
ical phenomenon at various scales. Most models of fluid 
systems rely on the Navier-stokes equation, which de- 
scribes the microscopic dynamics of the system. While 
the solutions to differential equations have been vali- 
dated many times in comparison with experiments, it 
remains challenging to predict the large-scale structure 
of the fluid flow in a turbulent regime. Recently topolog- 
ical techniques have been used to model and character- 
ize systems without reference to the microscopic dynam- 
ics. The present study develops characterizations that 
should offer insightful comparisons with other systems 
characterized using similar measures. In particular the 
evolution to power-law scaling has been of considerable 
interest 0, [H, [27 1 . 

While developing numerical models, one is not always 
able to match the initial conditions or exact boundary in- 
fluences that occur in experiments. The techniques out- 
lined in this paper could also be used to quantify global 
statistical measures as a tool to compare experimental 
and numerical results. The measures described here pro- 
vide average properties of the state of the system that 
complement other metrics currently in use [281 ] . Compar- 
isons based on topological measures have the advantage 
of being scale-independent. 

In spiral defect chaos, the extreme confinement of the 
fluid between the top and bottom plates primarily deter- 
mines the wave-number exhibited by the system. This 
shows up as finite and bounded distribution of power in 
Fourier space, corresponding to the wavelength of the 
convective rolls. We find that while the wave-vectors are 
constrained, the length scales of the rolls and topologi- 
cal measures characterizing the nearest-neighbor distri- 
butions have power-law distributions. 

The author would like to thank M. Gameiro, K. Mis- 
chaikow, N. Przulj and J. Scheel for useful discussions. In 
addition, the author is grateful to M. F. Schatz in whose 
laboratory the experiments were carried out. 
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